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5 ! Abstract 

Observed cusps with density profiles p oc r^^ or shallower, in the central regions 
of galaxies, cannot be reproduced in the standard Cold Dark Matter (CDM) picture 
l/-) ' of hierarchical clustering. Previous claims to the contrary were based on simulations 

■ with relatively few particles, and substantial softening. We present simulations with 
particle numbers an order of magnitude higher, and essentially no softening, and show 
that typical central density profiles are clearly steeper than p cx r~^. The observed 

■ shallower profiles may have formed through the smoothing effect of the spiral-in of 
central black holes in previous merger phases. In addition, we confirm the presence 
of a temperature inversion in the inner 5 kpc of massive galactic halos, and illustrate 
its formation as a natural result of the merging of unequal progenitors. 
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■ 1 Introduction 
> ' 

^ I Recent high resolution ground-based and Hubble Space Telescope observations (Lauer 

I et al. 1995) have revealed that elliptical galaxies do not have a constant-density core and 
the densities continue to rise until the resolution limit. In faint ellipticals, the density profile 
at the central region increases roughly as p oc r~^, and in bright ellipticals it increases as 
p oc or shallower (Merritt and Friedman 1996). 

Several simulations demonstrated that the dark matter halos formed through hierar- 
chical clustering were not well approximated by isothermal spheres, and better fitted by a 
Hernquist model (Hernquist 1990) or similar models with p ~ 1/r at the center (Dubinski 
and Carlberg 1991; Navarro, Frenk, and White 1996a, 1996b). In these simulations, halos 
have a power law index which changes from around —1 to —3 or —4 as the radius increases. 
Recently, Moore (1994) and Flores and Primack (1995) have argued that the gently ris- 
ing p(r) curves of dwarf galaxies are inconsistent with a central cusp, and therefore 
challenge the hypothesis of CDM halos. 

However, whether the results of numerical simulations are physically valid or are nu- 
merical artifacts remains unclear for the following reasons. First, no physical mechanism 
has been presented so far for the formation of cusps found in simulations (White 1996). 
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Secondly, since the mass resolution in these simulation was rather low, the central structure 
may have been affected by two-body relaxation effects (Quinlan 1996, Steinmetz and White 
1996, Fukushigc and Makino 1996). Thirdly, the central structures of halos formed in these 
simulations are strongly affected by the potential softening used, and show rapid change 
in the power index of p(r) around a radius not much larger than the potential softening 
length. 

In this letter, we have been able to separate physically real effects from numerical 
artifacts, through A^-body simulations of hierarchical clustering with a resolution far higher 
than those in previous work. The number of particles we used {N = 786, 400) is more than 
10 times larger than those used in previous simulations {N ~ 10, 000 — 30, 000), and our 
softening is significantly smaller (§2). 

Our results show that dark halos formed through hierarchical clustering typically exhibit 
an inwardly decreasing temperature structure in their inner regions, with a density cusp 
shallower than p oc r~^, but steeper than p oc (§3). We illustrate the origin of this 
behavior through an idealized merger simulation of different types of galaxies (§4), which 
enables us to give a physical explanation of temperature inversion (§5). 

2 Method 

Initial conditions were constructed following Dubinski and Carlberg (1991). We assigned 
initial positions and velocities to particles in a spherical region with a radius of 2Mpc 
surrounding a density peak selected from a discrete realization of a standard CDM model 
{Ho — 50km/s/Mpc and Q, — 1). The peak was chosen from a reahzation of the density 
contrast in an 8Mpc box, using COBE normalization. The peaks were found by smoothing 
the density with a Gaussian filter of radius 0.75 Mpc. 

We followed the evolution of a density peak through direct A^-body simulation. We 
added the local Hubble flow and integrated the orbits directly in physical space, with 
N — 786, 400, an individual particle mass of 4.0 x lO^M©, and a Plummer softened potential 
with a length of £ = 0.14kpc. We started the simulation at z ~ 46. We did not include tidal 
effects from outside our 2Mpc sphere, since we are mainly interested in the core properties 
within lOkpc, where tidal effects from outside 2Mpc scale are negligible. 

We used a 4-th order Hermite integration scheme (Makino and Aarseth 1992) with 
individual (hierarchical) timestep algorithm (McMillan 1986, Makino 1991a). 

For the force calculation, we used the GRAPE-4 (Taiji et al. 1996), a special- purpose 
computer designed to accelerate A^-body simulations using a Hermite integrator and a hi- 
erarchical timestep algorithm. The total system consists of 1692 pipeline processor chips 
dedicated to gravity calculation and has a peak performance of 1 Tffops. The calculation 
with N = 786,400 took about 180 CPU hours (5.7 x 10*^ particle steps) using 3/4 of the 
total system. The sustained speed of computation was 406 Gflops. Since the force calcu- 
lations on GRAPE-4 are effectively of double-precision accuracy (Makino ct al. 1996), our 
simulations exhibit much higher accuracy both for force calculations and orbit integrations 
than previous simulations. 
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3 Hierarchical Clustering Simulations 



Figure 1 shows the particle distributions in our simulation at the redshift z =(a)8.7, 
(b)5.1.and (c)1.8. Figure 2 shows the density and temperature structure of the halo at 
z — 1.8, well after most of the mergings has already taken place (results at z = are 
similar, but these earlier stages shows even less effects of two-body relaxation). In figure 2 
we can see a clear temperature decrease toward the center within 5kpc. This non-isothermal 
structure produces a density cusp shallower than p oc r~^. In our simulation, the structure 
outside 1 kpc is not affected by two-body relaxation after the formation of a large single 
halo at 2; ~ 3; only the central region within 1 kpc shows some expansion because of 
two-body relaxation effects. 

The large potential softening employed in previous studies has produced spurious struc- 
tures, which we have been able to reproduce, using additional simulations with larger soft- 
ening, leading to a central temperature decrease within a region a few times larger than 
the softening radius. The potential softening tends to produce a flat core on the scale of 
the softening length. Thus, unless we use point-mass particles, we always see a tendency 
for the power index of the density profile to approach zero at a radius comparable to the 
softening radius. 

After most of the merging has taken place, subsequent two-body relaxation effects due 
to small N continue to lead to spurious changes in the central region. For example, a 
density cusp might evolve to a flat core through gravothermal expansion (Quinlan 1996, 
Fukushige and Makino 1996). In our simulation, the local two-body relaxation times at 1 
and 5 kpc are 1.2 x 10^° and 2.3 x 10^^ years, respectively. Therefore, only the structure 
within Ikpc is somewhat affected by two-body relaxation in our simulation. In simulations 
with N ~ 10, 000 such as reported by Dubinski and Carlberg (1991) and Navarro et al. 
(1996a, 1996b), the local relaxation time at 5kpc was ~ 3.0 x 10^ years, implying that the 
effect of two-body relaxation completely dominated their results at that distance. 

4 Idealized Merger Simulations 

In order to illustrate the formation of temperature inversion after merging, we performed 
an idealized simulation in which we merged two equal-mass spherical halos with different 
central densities, a Plummer model and a King model with central potential l^o — 9. The 
result is shown in figure 3. The central temperature inversion is striking. 

The encounter between the halos is head-on, starting from rest at an initial separation 
of i? = 5 (in standard units with M = G = —4E = 1 where G is the gravitational 
constant, and M and E are the total mass and the initial total energy of a halo; cf. Heggie 
and Mathieu 1986). The total number of particles is N — 262,144 and the softening 
length is 1/256. We used Barnes-Hut tree code on GRAPE-4 (Makino 1991b) for the force 
calculations. The timestep is shared and constant {At = 1/512). 

In hierarchical clustering, temperature inversion takes place in the inner halos in a 
similar way. Figure 4 shows the evolution of temperature structure for the simulation of 
halo formation presented in section 3. We can see that the temperature of the outer halo 



3 



regions increases faster than the temperature of the central halo regions. 
5 Discussion 

In this letter, we have shown that dark halos formed through hierarchical clustering 
have a central density cusp shallower than p oc r^^, and a velocity dispersion that has a 
local minimum in the center, peaks at a distance of 5 ~ 10 kpc from the center, and then 
drops again in the outer halo. We offer the following interpretations, for the temperature 
and the density structures found. 

The occurrence of the striking temperature inversion can be understood as follows. 
In bottom-up structure formation, as in CDM hierarchical clustering, a typical halo is 
formed through repeated merging of smaller subclumps. Each time this happens, the 
less concentrated clump tends to be disrupted by the tidal field of the more centrally 
concentrated clump. The dense core of the latter survives the merging process more or 
less intact, and settles down at the center of merger remnant, with locally unchanged 
temperature. In contrast, the temperature of the merger remnant as a whole increases 
since the specific binding energy of the merger remnant is almost always larger than that 
of its progenitors. 

As a result, present-day halos tend to carry some memory of the temperature of their 
densest progenitor clump, which has set the temperature scale in the inner, and densest, 
regions. Since the bulk of the halo is affected more by subsequent mixing and heating, a 
central temperature inversion is naturally created. This formation process is evident, not 
only in our hierarchical clustering simulations, but also in our idealized simulation of a 
merger between two different halos, one with dense inner region, and another one with a 
much flatter core. 

The occurrence of steep density profiles, significantly steeper than p oc r~^, naturally 
follows from the same physical picture. According to our detailed A^-body calculations, in 
a purely CDM scenario, each final halo forms a repository for the small inner cores of the 
subclumps that have made up the final dark matter aggregate. 

Comparing our results with observations, we conclude that the shallow cusps of large 
ellipticals cannot have formed through the dissipationless processes that we have modeled 
here. Instead, we interpret the presence of density profiles shallower than p oc as strong 
evidence for the existence of central massive black holes, through the following mechanism. 
As an aftereffect of the merger of two black-hole containing galaxies, the spiral-in of those 
two massive black holes tends to smear out the central cusp that would have otherwise 
formed, as shown in detailed calculations of hierarchical merging of galaxies by Makino 
and Ebisuzaki (1996). 
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Figure 1: Snapshots of A^-body simulation of hierarchical clustering at redshifts (a) z — 8.7, 
(b)5.1, and (c) 1.8. The unit of length is 1 kilo parsec. 
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Figure 2: Density and temperature structures of the halo at 2; = 1.8 in the A^-body sim- 
ulation shown in Figure 1. The position of the center of the halo was determined using 
potential minimum and averaged physical values over each spherical shell. 
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Figure 3: The density and temperature structures of the merger remnant formed by the 
merging of two equal-mass clusters with different central densities. 
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Figure 4: Evolution of temperature structures of the halo at z = 8.7(triangle), 5.1 (square), 
and 1.8 (star) obtained by A^-body simulation illustrated in Figure 1. We made the tem- 
perature profile in the same way as in figure 2. 



9 



